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Abstract 

If one surrounds a black hole with a perfectly reflecting shell and adia- 
batically squeezes the shell inward, one can increase the black hole area A to 
exceed four times the total entropy S, which stays fixed during the process. 
A can be made to exceed 4S by a factor of order unity before the one enters 
the Planck regime where the semiclassical approximation breaks down. One 
interpretation is that the black hole entropy resides in its thermal atmosphere, 
and the shell restricts the atmosphere so that its entropy is less than A/ A. 
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The Generalized Second Law (GSL) of black hole thermodynamics states 
that the total entropy S does not decrease, and it further states that for Einsteinian 
gravity (to which this paper will be restricted, though the generalization to various 
other theories should be straightforward), 5* is the GSL entropy 

^GSL = ^^ + ^m, (1) 

where A is the total event horizon area of all black holes and S^^ is the entropy 
of matter outside the black holes. (I am using Planck units in which h, c, 4'7reo, 
Boltzmann's constant k, and the renormalized Newtonian gravitational constant G 
are all set equal to unity.) Although the Generalized Second Law has only been 
proved under restricted conditions, such as for quasistationary semiclassical black 
holes [0, it is believed to have greater generality, such as to rapidly evolving black 
holes. 

On the other hand, it is obvious that the GSL would not apply in the form above 
if the matter entropy were taken to be its von Neumann fine-grained entropy 

•S'vN = -^rplnp, (2) 

and if information is not fundamentally lost when a black hole forms and evaporates, 
since the initial and final fine-grained entropy of the matter could be zero when A = 
0, whereas the GSL entropy S'gsl = A/4: + S-^, which ignores quantum correlations 
or entanglements between the black hole and the matter outside, would be positive 
when (and only when) the black hole exists and has positive area. 

Nevertheless, under some suitable assumptions of coarse graining (such as ignor- 
ing entanglements between black holes and matter outside, and ignoring complex 
entanglements between different parts of the matter emitted from a black hole if the 
emission process is indeed described by a quantum unitary process), it has generally 
been assumed that it is a good approximation to take the total entropy S to be 
S'gsl, the sum of ^4/4 and S^- 

An implicit further assumption that is often made is that the matter entropy 
5*111 cannot be negative (as indeed it could not, for example, if it were given by 
= —trplnp). This assumption, plus the GSL, leads to the conclusion that the 
total entropy is bounded below by one-quarter the total event horizon area: 

\a < S. (3) 

Here I shall show that the inequality (^ can be violated. This violation can be 
interpreted as either a violation of the Generalized Second Law (if is assumed to 
be restricted to nonnegative values) or as an indication that the matter entropy 
must be allowed to take negative values in order to conform to the GSL. 
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Briefly, a violation of the inequality can be produced as follows: Take a 
Schwarzschild black hole of initial mass Mj and radius 2Mj ^ 1 (in the Planck units 
used herein) with negligible matter outside (e.g., before the black hole has had time 
to radiate significantly). Assuming the GSL for this initial state, the initial entropy 
Si is roughly Ai/4 = AirMf ^ 1, one-quarter the initial area of the hole, since the 
initial matter entropy is negligible in comparison. (If one considers as matter the 
thermal atmosphere that forms when the horizon forms in the near-horizon region 
r — 2M -C 2M, either the entropy of this atmosphere should be considered negligible 
if it is considered to be part of S^, or it should be considered to be part of the black 
hole entropy A/ 4; one gets too large a value for the total entropy if one counts both 
A/A and a large entropy associated with the near-horizon thermal atmosphere.) 

Now surround the black hole by a spherical perfectly reflecting shell at a radius 
that is a few times the Schwarzschild radius 2Mj of the black hole. This region will 
soon fill up with thermal Hawking radiation to reach an equilibrium state of fixed 
energy Mj inside the shell, but for Mj ^ 1, all but a negligible fraction (~ ^/Mf 
in Planck units) of the energy will remain in the hole, which can thus be taken still 
to have mass Mj. Outside the shell, one will have essentially the Boulware vacuum 
state with zero entropy (plus whatever apparatus that one will use to squeeze the 
shell in the next step, but this will all be assumed to be in a pure state with zero 
entropy). 

Next, squeeze the shell inward. If this is done sufficiently slowly, this should be 
an adiabatic process, keeping the total entropy fixed. Also, the outside itself should 
remain in a zero-entropy pure state, since the perfectly refiecting shell isolates the 
region outside from the region inside with its black hole and thermal radiation, 
except for the effects of the gravitational field, which will be assumed to produce 
negligible quantum correlations between the inside and the outside of the shell (as 
one would indeed get in a semiclassical approximation in which the geometry is 
given by a spherically symmetric classical metric). Some of the thermal Hawking 
radiation will thus be forced into the black hole, increasing its area. 

So long as the shell is not taken into the near-horizon region r — 2M -C 2M, 
the radiation forced into the black hole will have negligible energy and so will not 
increase the black hole area significantly above its initial value Ai. (Indeed, some 
of this tiny increase in the area just compensates for the tiny decrease in the black 
hole area when it filled the region r < rj with thermal radiation.) 

However, nothing in principle prevents one from squeezing the shell into the near- 
horizon region, where a significant amount of the near-horizon thermal radiation 
can be forced into the hole, increasing its mass M and area A = AnM"^ significantly. 
Since the entropy 5* should not change by this adiabatic process, it remains very 
nearly at Ai/A. Therefore, one ends up with a squeezed black hole configuration 
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with A > 4S ~ Ai, or total entropy significantly less than A/A. 

If the squeezing of the shell is accomplished purely by tensile forces within the 
shell that reduce its area, there is a minimum to the area of the shell (where the 
tensile forces needed to hold the shell in place diverge) and a maximum to the 
area of the squeezed black hole, such that A — AS cannot rise above some constant 
(depending on the radiation constant and hence on the number of matter fields 
present) times Mj (or some related constant times ^/S). However, if the squeezing 
is accomplished by applying radial forces to the shell, it can be squeezed past its 
minimum area to a sequence of configurations in which both its area and the area of 
the black hole inside increase by a fraction of order unity, with now the limitation 
being the onset of Planck-scale curvatures. To borrow the language from another 
field, when it is time to push, the black hole dilates. 

Perhaps the simplest way to incorporate these S < A/A configurations into black 
hole thermodynamics is modify the Generalized Second Law to state that 

'S'gSL = 'S'bh + Sm (4) 

does not decrease for a suitably coarse-grained nonnegative 5*111 and for a suitable 
definition of S'bh that reduces to ^4/4 (in Einstein gravity) when there are no con- 
straints on the near- horizon thermal atmosphere but which is less than A/ A when 
the atmosphere is constrained. One might interpret S'bh as arising entirely from the 
near-horizon thermal atmosphere, so that if the atmosphere is unconstrained in the 
vertical direction, its entropy is at least approximately A/ A. (There is no fundamen- 
tal difficulty in allowing that in this unconstrained case, S'bh might also have other 
smaller correction terms, such as a logarithm of the number of fields or a logarithm 
of A or of some other black hole parameter. It is just that in the unconstrained 
case, the leading term of Sbh should be proportional to A, and the coefficient should 
be 1/4, at least in Einstein gravity.) But if the near-horizon thermal atmosphere is 
constrained, it has less entropy. 

An alternative way to incorporate these S < A/A configurations is to retain the 
Generalized Second Law in the original form of Eq. (|I]), which is the special case 
of Eq. (^ in which Sbh = ^/4, but now to allow Sm to become negative when one 
squeezes the black hole. For example, one might use Eq. (|1|) not to define S'gSL 
in terms of A/ A and S'm, but instead to define S'm as the total entropy S'gSL minus 
the black hole entropy A/ A. (Of course, this procedure would make the GSL useless 
for telling what the total entropy is, so then Sgsl would have to be found by some 
other procedure.) 

An analogue in which such a definition would give negative entropies for some 
subsystems would be the case in which one used the von Neumann fine-grained 
entropy (|^) for the entropy of a total system (analogous to Sgsl) and for the entropy 
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of one subsystem (analogous to S'bh = A/ 4) and then simply defined the entropy 
of the second subsystem (analogous to Sm) to be the total system entropy minus 
the first subsystem entropy. If the two subsystems making up the total system 
have sufficient quantum correlations or entanglements, the entropy of the second 
subsystem thus defined can be negative. For example, if the complete system is 
in a pure state that entangles the two subsystems, the von Neumann fine-grained 
entropy of the total system would be zero, whereas the first subsystem would be in 
a mixed state with positive von Neumann entropy, so by the definition above, the 
entropy of the second subsystem would be the negative of that positive quantity. 

In pursuing this analogue, one could certainly take S'gsl to be the von Neumann 
fine-grained entropy of the total universe, though then under unitary evolution (e.g., 
no information loss in black hole formation and evaporation), this entropy would 
remain constant, so the Generalized Second Law would be rather trivial. (Or one 
could say the nontriviality is all in the fact that the evolution is unitary.) But 
the analogue would almost certainly break down for the black hole entropy S'bh if 
it is assumed to be A/ 4, since the von Neumann fine-grained entropy of a black 
hole would not in general equal A/ 4. One would expect it to be approximately 
A/ 4 when the black hole is maximally mixed for its area, but, at least under the 
assumption of unitary evolution, the actual von Neumann entropy of a black hole 
subsystem formed from some system of significantly smaller entropy (e.g., from the 
collapse of a star) would be expected to be much smaller than this maximum value 
of approximately A/ 4 until the black hole has emitted radiation (with which it thus 
becomes entangled) with entropy at least as great as A/ 4 of the remaining black 
hole @. 

For example, suppose we take a star of ten solar masses, or M ~ 10'^^ in Planck 
units. Its initial entropy will be of the order of the number of nucleons (of mass 
m„ ~ 10~^^) that it has. Si ~ M/nin ~ 10^^. If this star collapses into a black hole 
of the same mass (ignoring the mass ejection that would realistically take place), 
A/4 = 47rM^ ~ 10''^, but the von Neumann entropy would remain near 10^^, 21 
orders of magnitude smaller, for a very long time, ~ lO^MSi ~ 10^°° ~ lO^^yr, 
before the black hole increases its von Neumann entropy significantly above Si by 
emitting, and become entangled with, radiation of significantly more entropy than 
Si- If one squeezed the black hole so that A — 4S ^ S and assumed that the 
original GSL Eq. (|l]) were valid, one would need 5*111 ~ —10''^, which is not only 
negative but also is about 21 orders of magnitude in size larger than the relevant 
von Neumann entropies of the complete system and of the black hole that are of 
the order of Si. Therefore, it is rather hard to interpret such an enormous negative 

as arising from differences of these two von Neumann entropies, though it still 
might be possible in terms of field correlations across a region of the height of the 
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constrained near-horizon thermal atmosphere. 

Let us now try to estimate what the total entropy is of an uncharged, nonrotating 
black hole configuration of mass M and area A — IGttM^, in equihbrium with 
Hawking thermal radiation inside a perfectly reflecting pure-state shell of radius R 
and local mass /i, outside of which one has vacuum. We shall take a semiclassical 
approximation with a certain set of matter fields, which for simplicity will all be 
assumed to be massless free conformally coupled fields. Given the field content of 
the theory, the three parameters (M, R, fi) determine the configuration, though the 
entropy should depend only on M and R, since the shell and the vacuum outside 
have zero entropy. 

It is convenient to replace the shell radius R with the classically dimensionless 
parameter 

W ^ f . (5) 

which would be if the shell were at infinite radius (though before one reached this 
limit the black hole inside the shell would become unstable to evaporating away) 
and 1 if the shell were at the black hole horizon (though in this limit the forces on 
the shell would have to be infinite). Then we would like to find S{M, W). 

If W is neither too close to nor to 1, the entropy will be dominated by A/4 = 
47rM^. The dominant relative correction to this will come from effects of the thermal 
radiation and vacuum polarization around the hole and so would have a factor of h if 
I were Tising gravitational units (c — G = 1) instead of Planck units {h — c — G = 1). 
In gravitational units, h is the square of the Planck mass, so to get a dimensionless 
quantity from that, one must divide by (or by R^, which is just AM'^/W'^ with 
W being of order unity); for the free massless fields under consideration, there are no 
other mass scales in the problem other than the Planck mass. Therefore, in Planck 
units, the first relative correction to A/ A will have a factor of and hence give 

an additive correction term to AnM"^ that is of the zeroth power of M, a function 
purely of W. One might expect that if one proceeded further in this way, one would 
find that the entropy S is given by AnM'^ times a whole power series in with 
each term but the zcroth-order one having a coefficient that is a function of W. If 
we had been considering the possibility of massive fields, then these coefficients of 
the various powers of would not be purely functions of W but would also be 

functions of the masses of the fields. However, for simplicity we shall consider only 
the free massless field case here. 

In fact, I shall consider only the first two terms in this power series, 

S(M,W) = A7rM^ + fi{W) + f2{W)M-^ + fs{W)M-^ + --- 

^ A7rM^ + fi{W). (6) 
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The function fi{W) will depend on the massless matter fields present in the theory, 
most predominantly through the radiation constant 



«r- = ^K + ^n/), (7) 

where nj, is the number of bosonic helicity states and is the number of fermionic 
helicity states for each momentum. It also proves convenient to define 

^ Qr ^ rib + jrif 

~ 3847r3 11 5207r ' ^ ^ 

which makes the entropy density of the thermal Hawking radiation far from the hole 
(when R:$> M or W <^1) simply a/M^, and also to set 

MW) = -327ras{W), (9) 

where s{W) depends (weakly) only on the ratios of the numbers of particles of 
different spins and so stays fixed if one doubles the number of each kind of species. 
Then my truncated power series expression for the entropy of an uncharged spherical 
black hole of mass M at the horizon (and hence horizon radius 2M and horizon area 
AttM"^) surrounded by a perfectly refiecting shell of radius R — 2M/W is 

S{M, W) « AttM^ - 32TTas{W) = \a - 32TTas{W), (10) 

Now I shall evaluate an approximate expression for s{W) when the perturbation 
to the Schwarzschild geometry is small from the thermal radiation inside the shell 
and from the vacuum polarization inside and outside the shell. There will be an 

additive constant to s{W), giving an additive constant to the entropy, that I shall 
not be able to evaluate, but for simplicity and concreteness I shall assume that 
s(l/2) = 0, so that the entropy is A/A when the shell is at W = 1/2 or i? = 4M. 

First, I shall ignore the Casimir energy and related effects of the shell itself 
on the fields. I would expect that these effects would give additive corrections to 
s{W) that are of order W or smaller (and so never large compared with unity), 
whereas the leading term in the perturbative approximation for s{W) will go as 
(proportional to the volume inside the shell) ior W <^ 1 (shell radius R » 
2M) and as 1/(1 — W) (inversely proportional to the redshift factor to infinity) for 
1 — W<^1 (shell radius relatively near the horizon), so one or other of these leading 
terms will dominate when W is near or 1. Therefore, I shall take the stress-energy 
tensor inside the shell to be approximately that of the Hartle-Hawking state in the 
Schwarzschild geometry, and that outside the shell to be approximately that of the 
Boulware vacuum. 
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The first part of tlie analysis will be done in a coordinate system {x^, x^, x"^, x^) = 
{t, r, 6, 0) in which the spherically symmetric classical metric has, at each stage of 
the process, the approximately static form 



ds'' = -e^'f'dt^ + U-^dr^ + r\de^ + sin^ Odip^) (11) 



with 



and 



with 



e'f = e'^U (12) 
2m 

U=l = l-w 13 

r 

yj^fUl = l-U =1- (Vr)2. (14) 
r 

Here 0, U, ip, ^i and w are all functions of the x^ = r coordinate alone, although 
they also have a global dependence on the black hole mass M (the value of r/2 
where e^'^ = 0), and (for r > R) on the radius r = R = 2M/W of the shell and on 
the total rest mass energy ^ of the shell. The Einstein equations then give 

^ = ATTr(p + P)U-^ (15) 
dr 

and 

dm , n , . 

— = Anr^p, 16 
dr 

where 

P = -(To°(r)) (17) 
is the expectation value of the energy density in the appropriate quantum state, and 

P=(T^(r)) (18) 

is the corresponding expectation value of the radial pressure, both functions of r. 
The functional form of the expectation value of the tangential pressure {T2{r)) = 
(T|(r)) would then follow from the conservation of (T^^) but will not be explicitly 
needed in this paper. 

Since we are assuming that the state of the quantum fields inside the shell (r < R) 
is the Hartle-Hawking [^] thermal state, for r < R we have 

p = pH{M,r)^^^^eH{w) (19) 

and 

ct 

P = PH{M,r) ^ ^^-^^pj,{w), (20) 
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where on the extreme right hand side of each of these two equations I have factored 
out the dependence on the black hole mass M from that on the classically dimen- 
sionless radial function w = 2m/r, thereby defining two classically dimensionless 
functions of w, Suiw) and ph{w). 

Similarly, we are assuming that the state of the quantum fields outside the shell 
(r > R) is the Boulware [|5| vacuum state, so ioi r > R we have 

and 

P = PB{M^,r) ^ ^^^p^{w), (22) 

thereby defining two new classically dimensionless functions of w, SBiw) and pb{w). 
Here 

Moo = m{r = oo) (23) 

is the ADM mass at radial infinity. 

In some cases below we shall assume that there is some extra apparatus in the 
region r > R holding the shell in. If so, its energy density and radial pressure can 
simply be included in ps and Pg- In any case, we shall assume that whatever is 
outside the shell is in a pure state with zero entropy. 

Below we shall also need the vacuum polarization part of the stress-energy tensor 
inside the shell, whose components I shall denote by 

pv{M, r) = ph{M, r) - pt(M, r) 

3a , , 3a , , , 

^viw) = 7;^r-TT7ieHiw) - Sriw)) (24) 



32nM^ ' ' 327rM4 
and 

Pv{M, r) = Ph{M, r) - Pt{M, r) 

= rT7Py('U^) = tt7[Ph(w) — pt(w)), (25) 

327rM4^^ ^ 327rM^^^^ ; v \ j 

where pT and Pt denote the components of the thermal parts. 

I shall assume that the vacuum polarization part is what the Boulware state 
would give if one had it inside the shell, so that pv and Py have the same dependence 
on the local mass m(r) and radius r as ps and Pb do outside the shell (when there 
is no extra apparatus there). In the first-order (in a/M"^) perturbative calculation 
being done here, the expectation value of the stress tensor is already first order 
(except possibly for that of the shell), so its functional dependence on m can be 
replaced by its dependence on its zeroth approximation, which is the black hole 
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mass M for r < R and the ADM mass Moo for r > R. Therefore, to sufficient 
accuracy for our purposes, pv and Py can be evaluated by using Eqs. ( pT]) and 
([2^ ) for pb and Pb with the ADM mass Moo, which is approximately the value of 
the local mass m(r) anywhere outside the massive shell, replaced by the black hole 
mass M, which is approximately the value of m(r) anywhere inside the shell. In 
particular, this implies that we can use 

sviw) = €b{w) (26) 

and 

py(w)=pBiw) (27) 

For explicit approximate calculations, it is useful to have explicit approximate 
formulas for these various components of the stress-energy tensor, given by the 
equations above from the six functions Suiw), SBiw), Sxiw) = euiw) — eB^w), 
Ph{w), pb{w), and Pt{w) = puiw) —pb{w). For simplicity and concreteness, I shall 
use those obtained for a conformally invariant massless scalar field in the gaussian 
approximation which gives 



, , 327rM^ (SttMY 1-U-3wYw^ 
eniw) = — ph = ^ ^-Ph^ ^ 2Aw 

l + 2w + 3w^ + Aw^ + 5w^ + 6w^ -33w^, (28) 



and 



, , 327rM^ (SttMoo)^ -U-3wfw^ ^ . 

^ l + 2w + 3w'^ + 4:w^ + 5w^ + 6w^^ - 33w^ (29) 



1 — w 



|2 



eT{w)=eH{w)-eB{w)^-^j-^-^ = ^, (30) 



, , 327rM4 ^ (SvrM)^ ^ 1 - (4 - 3wYw'' , 
Ph{w) = Ph = \ ' Ph ~ r^ \i — + 24w^ 

l + 2w + 3w^ + Aw^ + + 6w^ + 15w^ (31) 



a Sttr [1 — wY 

^ +l + 2w + 3w'^ + Aw^ + bw^ + Qw^ + 15u7^ (32) 



2 



{l-w) 



Pt{w) = Ph{w) - Pb{w) ^ = (33) 
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Note that this approximation gives 



PT ^ Wt ^ a,Ti*_i, (34) 

just hke thermal radiation in flat spacetime, where Tiocai is the local value of the 
Hawking temperature, 

T,ocai^^(l-— )-^/^ (35) 

Now we use the Einstein equations (|T5|) and (|16D with the appropriate p and P on 
the right hand side, and with the metric function U there taking on its approximate 
Schwarzschild form, 1 — 2M/r for r < R and 1 — 2Moo/r for r > R. 

We also need to consider the effect of the shell, which has a surface stress-energy 
tensor with components 

and 

Si = Sl = — ^, (37) 

where p is the total local mass of the shell, the shell area AttR'^ multiplied by the 
local mass-energy per area —Sq as seen by a local observer fixed on the shell, and F 
is the local total tensile force pulling together the two hemispheres of the shell, the 
circumference 27rR multiplied by the local surface tension (tensile force per length) 

02 _ c3 
^2 ~ *-'3- 

If one integrates the Einstein equations (|15D and (|16]) through the shell and uses 
the conservation law for the stress-energy tensor, one get the junction conditions ^ 
in the static case that 

= i?(f/y2 _ f/1/2) (3g) 

and 

8F = + (1 + 87ri?2p_)t/rV2 _ (^ ^ 87ri?2p^)t/-V2^ (39) 
R 

where 

U-^l-'-f (40) 

is the value of U just inside the shell (r = R—), where the local mass function m 
takes on the value M_, and 

U.-l-'-f (41) 

is the value of U just outside the shell (r = R+), where the local mass function m 
takes on the value M+. Similarly, P_ and P+ are the expectation values of the radial 
pressure of the respective quantum states just inside and just outside the shell. 

Thus we have at least five relevant masses for the configuration: the black hole 
mass M = m{r = 2M), the mass M_ = m(r = R—) just inside the shell, the 
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local mass (or local energy) fi of the shell itself at radius r = R, the mass M+ = 
m(r = R+) just outside the shell, and the ADM mass M^o = rnir = oo) at radial 
infinity. Since the stress-energy tensor inside the shell is that of the Hartle-Hawking 
state determined by M and r, M_ is a function of M and R. Similarly, since the 
stress-energy tensor outside the shell is that of the Boulware state determined by 
Moo and R (at least when we do not have an extra apparatus there to hold the shell 
in place), M+ is a function of Moo and R. The junction condition ( ^Sf ) then gives /i 
as a function of M_, M+, and R, and hence as a function of M, Moo, and R. One 
can in principle invert this to get Moo (and hence all the other masses as well) as a 
function of M, R, and n, or to get M and all other masses as a function of Moo, R, 
and /i. The main point is that if we just have a black hole with the Hartle-Hawking 
thermal state inside a shell, and the Boulware vacuum state outside the shell, the 
semiclassical configuration (for fixed field content of the quantum field theory) is 
determined by three parameters, though only two of them (say M and either R or 
W = 2M/R) are relevant for the entropy which resides purely inside the shell. 

To evaluate the function s(W) in the truncated entropy formula (|TO|), I shall 
consider an adiabatic process of slowly squeezing the shell, keeping the total entropy 
constant and thereby getting 

ds _ M dM 

dW~A^dW ^ ' 

during this process. Since this process is not strictly static, one cannot use precisely 
the static metric (|l]) with and U (or ip and m) that are purely functions of r 
and obey the static Einstein equations (|^) and (p!6D. However, one can consider 
a quasi-static metric in which and U (or ip and m) have a very slow dependence 
on the time coordinate t and the Einstein equations are only slightly different from 
Eqs. (0) and (|1§. 

The specific calculation which I shall do will be to have the shell squeeze itself 
inward by using its own internal energy, so that no apparatus is used outside the 
shell to push it inward, and that region has only the Boulware vacuum polarization. 
The contraction of the shell is assumed to be so slow that it does not excite the 
vacuum outside it but rather leaves it in the Boulware vacuum state with constant 
Moo. However, as the shell moves in, it is enlarging the Boulware state region, so 
effectively the shell must be creating a larger volume of vacuum with its vacuum 
polarization. This means that in the slowly inmoving frame of the shell, there is 
a flux of energy from the shell into the Boulware region, needed to enlarge the 
Boulware region while keeping it static where it already exists. [For the stress- 
energy tensor components of the Boulware vacuum given by Eqs. (^9]) and (0), 
this energy influx into the Boulware region is actually negative, so it increases the 
energy of the shell as it moves inward.] 
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Similarly, if the inside of the shell were also vacuum that did not get excited by 
the adiabatic contraction of the shell, there would be a swallowing up of part of the 
vacuum region by the shell as it moves inward. This would give a flux of (negative) 
energy from the vacuum inside into the shell, decreasing its energy. Surely this flux 
into the shell also exists even if the inside is not vacuum, and I assume that it is 
given by the vacuum polarization part of the actual stress-tensor there, which I take 
to be approximately that of a Boulware state with the same m and r. The remaining 
part of the total stress-energy tensor there, which I am calling the thermal part, and 
which is given approximately by Eq. (p4D, should simply be reflected by the shell 
and not give an energy flux into it (in the frame of the slowly contracting shell), 
though it will contribute to the force that needs to be counterbalanced very nearly 
precisely to obey the static junction equation ( p9D to high accuracy in order that 
the shell not have any significant acceleration relative to a static frame. 

In other words, I am assuming that if a shell moves inward through a static 
geometry, the vacuum polarization part of the stress-energy tensor will stay static, 
with T° = -pv{M,r), = Pv{M,r), and = T° = inside the shell, and with 
TO = -pB(Moo,r), = PB{M^,r), and = T° = outside the shell. Then as 
the shell moves through this static stress-energy tensor, in the frame of the shell, 
there will be fluxes of energy into or out from the shell on its two sides. In constrast, 
I am assuming that the thermal radiation part of the stress-energy tensor will be 
perfectly reflected by the shell, so that in the frame of the shell it will give no energy 
fluxes into or out from the shell. 

There is a modiflcation of this picture that occurs when the inward motion of the 
shell squeezes thermal radiation into the black hole so that its mass goes up. While 
the hole mass is increasing, the vacuum polarization inside the shell is not quite static 
but instead has small Tq and T° terms that, for sufficiently slow adiabatic processes, 
are proportional to M, the coordinate time derivative of the black hole mass M. In 
the present calculation, in which the shell is squeezing inself inward by using its own 
internal energy, the ADM mass Moo stays fixed, and so the vacuum stress-energy 
tensor outside the shell stays static during the process, under my approximation of 
neglecting Casimir-type boundary effects of the shell itself on the quantum field. 
For a sufficiently slow inward squeezing of the shell, the Tq and terms inside are 
small, but over the correspondingly long time of the squeezing they contribute an 
effect on the energy balance of the shell that is not completely negligible when one 
contemplates squeezing the shell to a final position very near the black hole horizon. 
(My original neglect of these terms caused me considerable confusion during early 
stages of this work and lectures I gave about it.) 

My present procedure for calculating the small and terms inside the shell 
is to assume that the shell squeezing, and all consequent processes, occur so slowly 
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that Tq and are given to high accuracy by the same functions of M and r as they 
are when the geometry is static, namely —pv{M,r) and Pv{M,r). Then I assume 
that the vacuum polarization part of the stress-energy tensor is itself conserved away 
from the shell, so one can use the conservation of its energy to deduce the radial 
derivative of e^r'^T^. 

In particular, if we let the vacuum polarization part of the stress-energy tensor 
have the component 

, _ aMw' 

^°"47rM4^ ^' ^^^^ 
with the factors chosen so that / is a function purely of w, then Tq.^ = becomes 

For the region inside the shell with r not too much larger than 2M, one has m ^ M 
and e'^ ~ 1 (possibly after suitably normalizing the time coordinate t) . Then if one 
uses Eqs. (p^)-(|27|), one can rewrite Eq. (|4^ ) as 

(45) 



df 3w d f6B\ <i£B+PB 



dr 4 dw Kw'^J 8w'^(l 



w 



Given the functions esiw) and p_b(w), e.g., as given by Eqs. (^) and (^) from the 
gaussian approximation for a conformally invariant massless scalar field, one can 
integrate Eq. ( ^5]) to obtain f{w) up to a constant of integration. Although the 



constant of integration is not important, it can also be fixed by assuming that an 
observer that remains at fixed w = 2m /r as m changes sees in its frame no energy 
flux in the limit that w is taken to unity, which implies that the flux of vacuum 
polarization energy through the horizon is taken to be zero. 

After one calculates the vacuum polarization part of the stress-energy tensor, 
which gives — Tq = Pb + Pb and = outside the shell and — Tq=pv + Pv 
and Tq as given by Eq. (|^) inside the shell, one can then calculate the fluxes of 
energy out from and into the shell and insert these into the conservation equations 
for the surface stress-energy tensor of the shell. For a very slowly expanding or 
contracting shell, one finds that 

dp = AFdR + A'KR'^dR[{pB + Pb)U+^''^ - {pv + Pv)UI^^^] - A-KR^T^UZ^'^dt. (46) 

The first term on the right hand side is the work done by the tensile force within 
the shell, and the remaining terms are the energy input from the vacuum stress- 
tensor components ps and Pb just outside the shell and the vacuum stress-tensor 
components pv, Pv, and Tq just inside the shell. 
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One now combines this local energy conservation equation for the shell with the 
static junction equations (|38|) and (^) that should still apply to high accuracy in 
this slowly evolving situation to keep the shell radius from accelerating too rapidly. 
When one also combines this with the integrals of Eq. (|16D, 



M_ = M + / 4:7rr^drpH, 

POD 

JR 



one finds 



1 - 



M2 



dM ^ -ATxR^dR{pT + Pt) 



(47) 
(48) 

(49) 



during the adiabatic contraction of the shell, which, up to the small correction factor 
involving /, is precisely what one would get in flat spacetime from adiabatically 
contracting a ball of thermal radiation. 

Next, we can use the fact that R = 2M/W to derive that 



dW 
dM 



M dR" 
'RdM 



2 
R 



1 + 



M(l -4a//M2 



(50) 



where / and pt + Pt are to be evaluated at r 
into Eq. (021) then gives 



AnR^ipT + Pt) 
R OT w ^ W. Inserting this back 



ds 
dW 



3£b + Pb 



4a 



3gi? + Vb 



f 



(51) 



For massless particles of any spin, it should be a good approximation to take pt + 
Pt ~ (4/3)0^71^^.3^1 in terms of the local temperature Tiocai, and this implies that 



+ Pb 



4/(1- 
ds 



Wf 



so 



dw w^{i - wy 



1 + 



4a 
M2 



w^{i-wy 



-f 



(52) 



If we omitted the / term from the radial flux of vacuum polarization energy 
when M changes, as I indeed first erroneously did, then the factor inside the curly 
brackets above would diverge as one approached the horizon, where W = 1. This 
implies that the reciprocal of this factor would cancel the divergence in the factor 
before it, so ds/dW would stay finite all the way down to W = 1, and one would 
find that the increase of one-quarter the area over the entropy would be limited to 
an amount of the order of aM. For large M this is large in absolute units, but it is 
always much smaller than the entropy itself, which is of the order of AttM"^. 

However, one can use the fact that the regularity of the Hartle-Hawking stress- 
energy tensor at the horizon imphes that pn + Ph, and hence Sen + Ph, must go 
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to zero at least as fast as 1 — w as one approaches the horizon. (This is easiest to 
see in the Euchdean section with imaginary time t, on which for fixed coordinates 9 
and ip, the horizon is at the center of a regular rotationally symmetric two-surface 
with angular coordinate int with k ^ 1/(4M) being the black hole surface gravity 
and with the radial distance being roughly 4M-\/r^--^ when 1 — w ^ 1. Then 
Ph = Tl is the pressure in the radial direction, and —pn = Tq is the Euchdean 
pressure in the Euclidean angular direction, and regularity at the origin demands 
that the difference go to zero at least as fast as the square of the radial distance 
from the origin.) Then one can show that / cancels the divergence in W~^{l — W)~^ 
so that — — / stays finite as one approaches the horizon. In fact, if 

one chooses the constant of integration of / so that the flux of vacuum polarization 
energy through the horizon is zero as M is slowly changed, then — W)^^ — / 

actually goes to zero linearly with 1 — W as one approaches the horizon. For example, 
using this constant of integration and the gaussian approximation for Ssb and pb 
leads to 

ds 1 ( 4a 



dw ^ w\i - wy 



|l + J^^C^ - + 3W + 6W^ + 2W^ + 7W* + 13W' 

(53) 

Therefore, we see that the correction term that is first order in a/M'^ inside the 
curly brackets of Eqs. (|5lD -(^3|) does not diverge as one takes 1 — W to zero but 
instead always remains small. Therefore, we can drop it (as we have also neglected 
other finite corrections that are linear in a/M"^) and integrate the zeroth-order part 



of Eq. (|53|) to get an explicit formula for s{W): 

_ dw _ 1 W 1 L_A ^ 

As discussed above, I arbitrarily chose the constant of integration of this integral to 
make s{W) = at W = 1 /2 or R = 4M, but this is not likely to be valid, and there 
are also Casimir energy effects from the shell and corrections to Eqs. ( |5(]| ) and ( PHj ) 
that would give correction terms at least of order W and likely also of the order of 
a constant and of order 1/W. 

Finally, we can insert this form for s{W) into Eq. ([To|) to get 

S{M,W) = ^A[l-^s{W)]=47iM^ - 327ias{W) 

- ^-^^-^2.«[^ + 41n^-^-^ + 0(^)] 

^ 4nM'-32na[ ^ + 41n^^^^^ ^ _ J^l, (55) 

^R-2M R-2M 24M3 AM^^' ^ ' 



where after the last approximate equality I have dropped the terms in Eq. (^J) that 
I suspect are always dominated by corrections to my approximations that I have not 
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included. Although I have retained four terms from s{W), only the first two terms 
should be kept when 1 — W = 1 — 2M/R -C 1 (shell very near the horizon), and 
only the next two terms should be retained when W = 2M/ i? -C 1 (shell very large 
compared with the black hole). 

The result indicated by Eq. (0) is precisely the same that one would obtain by 
taking the geometry to be Schwarzschild with a thermal bath of radiation with local 
Hawking temperature 

and entropy density {4:/3)arTj^c^i, and then taking the total entropy to be AnM"^ 
plus the entropy difference between that inside the shell at radius R and that inside 
the radius 4M. If one naively integrates this assumed entropy density all the way 
down to the horizon, one would get a divergence, but one can take the attitude that 
this divergence is regulated so that the entropy in this thermal atmosphere below 
some radius like AM (the precise value of which doesn't matter much, since the 
assumed entropy density is this region is of the order of a/M^) is the black hole 
entropy S'bh ~ A/4 = AirM"^. Then one can say that if the shell is put at a much 
larger radius, the entropy of the thermal Hawking radiation outside 4M or so would 
be matter entropy S^a that would add to 5'bh, which is certainly an uncontroversial 
assumption. 

What I have found from my consideration of having the shell squeezed in adi- 
abatically is that if the shell is put much nearer the horizon than a radius of AM 
or so is, then the entropy is correspondingly less than the usual black hole entropy 
S'bh ^ A/A = AnM'^. Because the thermal atmosphere is restricted from filling up 
the region to AM or so, it does not have the entropy needed to make the total 
entropy as large as A/ A. 

The next question is the range of W over which one would expect that Eq. ( ^5]) 
is approximately valid. For very small W or very large R, one essentially has a black 
hole of mass M surrounded by a much bigger volume, V ~ AttR^/3, of radiation in 
nearly flat spacetime with Hawking temperature l/(87rM)~^, energy density roughly 
3a/(327rM^), and entropy density roughly a/M^. The dominant term for the total 
energy of the radiation is Er ~ aR^ / {8M^) , and from Eq. (|55|), the dominant term 
for the total entropy of the radiation is Sr ~ AnaR^ / {3M^) . This agrees with the 
standard expression for the entropy of thermal radiation of energy Er in a volume 
V, 

Sr = liarVy/'E'r^' ~ fa'^^SRErf^ (57) 
For fixed total energy Moo = M + Er <ti the total entropy 

S « AttM^ + 5, ~ 47r(Moo - Erf + ^a^/\mErf'^ (58) 
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is indeed extremized for 

^'^8M^^ 8{M^-Er)^' ^^^^ 

but the extremum is a local entropy maximum if and only if 5Er < Moo or AEj. < M 
@, which implies that one needs R < (2M^ / aY^^ or 

for thermodynamic stability. 

For smaller values of W (larger values of R), the radiation energy E^- is more 
than 20% of the total available energy Mqo (assumed to be held fixed), and then 
if the black hole emits some extra radiation and shrinks, it heats up more than 
the radiation does, leading to an instability in which the black hole radiates away 
completely. On the other hand, if the black hole absorbs some extra radiation, it will 
grow and cool down more than the surrounding radiation, therefore cooling down 
more and absorbing more radiation, until the radiation energy drops to the lower 



positive root of Eq. (p9D, which is less than 0.2Moo and hence is at least locally 
stable. 

At the opposite extreme, the question is how small 1 — W can be. Here the 
fundamental limit is the Planck regime, which is the boundary of the semiclassical 
approximation being used in this paper. The Boulware vacuum energy density pb 
just outside a massless shell (so that the mass just outside, M+, is very nearly the 
same as the black hole mass M; for positive shell mass fi, pb would have an even 
greater magnitude) is, for very small U = 1 — W, ps ~ — 3a/(327rM^?7^). Suppose 
the semiclassical theory is valid until the orthonormal Einstein tensor component 
Gq = —SnpB ~ 3a/(4M^[/^) reaches a maximum value of, say. Cm, which would 
be expected to be of order unity (orthonormal curvature component of the order of 
the Planck value). This gives the restriction 

For U = 1 — W <^ 1, the spatial distance from the shell to the horizon is 
D ~ 4Mf/^/^, so this restriction on U gives a minimum distance the shell can be 
from the horizon: 

in Planck units, as all quantities are in this paper unless otherwise specified. 

If we combine the restriction (|6lD with the lower bound on W from Eq. (^) and 
re-express the combined restriction as a restriction on the radius R of the shell, we 
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get 

Alternatively, in terms of the distance D of the shell to the horizon (which is D ~ i? 
for R > 2M), we get 

If we now insert the restriction (|6TD or (|62D into the asymptotic form of the total 
entropy ( |55D for t/ = 1 — ^ 1, which is 

5(M,iy) ~ 4vrM2 - ^ ~ 4vrM^ (l - 1^ 

~ 47rM^ - ~ 47rM^ (^1 - j , (65) 

we get the limitation 



S{M,W) > MP ( 1 - 16^/^ 1 = AnM' | 1 - + ^^/) 1 • (66) 

This can be re-expressed as a limitation on how much the area A of a black hole 
can exceed four time the entropy, 45*: 



Therefore, unless we have N = rib + 7nf /8, the effective number of one-helicity 
particles, comparable to or greater than 1357r/Cj\/ ~ 424/Cm, the fractional increase 
of the black hole area A above 45" is restricted to be rather small, though even just 
N = 4 from two-helicity gravitons and photons would give a fractional increase of 
about 10% if the curvature limitation Cm is one in Planck units. 

However, this does raise the interesting theoretical question of what would hap- 
pen in a theory in which the effective number of particles is so large that 
A^C'Af/(1357r) is bigger than unity. Naively it would then appear that, without run- 
ning into excessive curvatures (i.e., Gq > Cm), one could put the shell sufficiently 
close to the black hole that the total entropy, given approximately by Eq. (|6^), 
would be negative, which is surely nonsense. Of course, one could never get to such 
a configuration by adiabatically compressing the shell, since one would then have 
started out with positive entropy that would not decrease (though it does appear 
that one could in principle push hard and far enough on the shell that the black 
hole could dilate to an arbitrarily large radius for the fixed initial entropy, which by 
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itself would be rather remarkable). However, one could imagine constructing such a 
shell in place and then evacuating the thermal radiation from above it, which would 
seem to leave behind a black hole configuration of negative entropy. In our world 
this possibility might be excluded by the limited number of particle species that 
exist and contribute to the thermal Hawking radiation, but one would like to see a 
direct argument of why negative total entropies of black hole configurations cannot 
be achieved even in a model universe in which one had a huge number of fields. 

The first guess that came to me off the top of my head is that of course no physical 
shell can be perfectly reflecting. A partially reflecting shell should be able to squeeze 
enough of the Hawking radiation into a black hole to raise its area above four times 
the entropy, but the transmission will put a limit on how effective this process can 
be. If there are more species of particles that can be partially transmitted, the 
shell may become less effective in increasing the black hole area, and conceivably 
this could offset the increase in the otherwise theoretically allowed fractional area 
increase from the increase of the number of species. However, this is just a wild 
guess, and so the problem will be left as an exercise for the reader. 

One refinement of the results above that should be explained is that although 
I have used a Schwarzschildean coordinate system (|ll|) in the analysis described so 
far, this system breaks down when one follows the Boulware vacuum sufficiently far 
inward. In particular, the Schwarzschildean radial coordinate r, which is l/(2-7r) 
times the circumference of the two-sphere, is only a good coordinate when it de- 
creases monotonically inward. This indeed occurs when U = 1 — 2m/r = {dr/dDy 
remains positive, where D is the proper radial distance. But because the Boulware 
vacuum energy density p_B is negative for sufficiently small U, as one moves inward 
with r initially decreasing, the mass function m(r) increases rather than decreases, 
and eventually one reaches a radius where 2m = r and hence U = 0. In this case, 
"fl'oo = e^*^ remains positive (which implies that e"^"^ = e^'^U~^ diverges, so ip is no 
longer a good metric function either), which means that one has not reached the 
horizon, but rather just a location where the gradient of the Schwarzschildean ra- 
dial coordinate vanishes. Inside this point, the gradient of r reverses sign, so now r 
increases as one moves inward. The Einstein equation (p!6| ) implies, since the energy 
density remains negative, that the mass function m now decreases as one moves 
inward. 

Since at small radii the Boulware vacuum polarization gives a stress-energy tensor 
that is dominated by contributions that look like thermal radiation but with the 
opposite sign (e.g., an isotropic pressure that is one third the energy density, both 
of which are negative), one can incorporate an approximation for its back reaction 
on the metric simply by solving the Tolman-Oppenheimer-Volkoff equations for such 
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a thermal fluid with 

If one defines the function 

ly = 1 + Snr'^PB (69) 

(not to be confused with the previous use of W as 2M/R), then one gets the two 
differential equations 

U{4-U -3W)dW = 2{W -2U){l-W)dU 

= 2{W -2U){1-W){A-U ~3W)dr/r. (70) 

One can match with the Schwarzschild metric, slightly perturbed by the Boul- 
ware vacuum polarization, at y/a/Moo <^ f/ ^ 1, where 1 — W ^ aM~^(l — f/)"^f/~^ 
and r ^ 2Moo(l — U)^^. Then as one integrates Eq. ([70| ) inward, initially r, U, and 
W all decrease, until U and W simultaneously go to zero. This is a singular point of 
Eq. (^^, but one can easily see that U goes to zero quadratically with W, so that 
as W crosses zero and becomes negative, U becomes positive again. If we define the 
new variable 

Y = ^= ^ " '^'^1'' m\ 

W l + 87rr2p' ^ ' 

we find that it decreases monotonically as we go inward, starting at the small positive 
value Y = Yq <^ 1 m order that one be in the regime where Eq. (^) is valid, but 
then with Y going negative where U goes to zero quadratically with W . 
In terms of Y and W , the first differential equation of ([70|) becomes 

dW_ ^ 2{l-2Y)W{l-W) 

dY Y[{l + 2Y){2-W) -WW)]' ^ ' 

Although this equation is also singular at F = and ly = 0, YdW/WdY = 1 there, 
so W just passes through zero linearly with Y, with a positive coefficient that one 
can calculate is roughly M^j ^fa. Then one can easily see that the right hand side 
of Eq. (^) is always positive, since 1 — TV > and 1 — 2Y > everywhere where 
Eq. ([68| ) is valid, so Y and W are both monotonically decreasing variables as one 
moves inward through the negative pressure (and negative energy-density) vacuum 
Boulware stress-energy tensor. 

Eq. (^) is not separable, but by pulling out the separable parts, one can get 



W 

oc exp 



(8 - 7W)dY 



{2 - W) + {A - 5W)Y 



(73) 



For the initial conditions above with y/a/Mao -C 1 one can now integrate Eq. 
approximately in the various regimes for W and match them: 
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Initially, W is near 1, and so 

W M, 



oo 



Y). (74) 



In particular, the right hand side is correct to first order in Y when i/a/Moo <^ Y. 
Then as W goes from near 1 to much less than —1, Y is so small and changes so 
little that the integral in Eq. ( [73| ) has only a negligible contribution. Therefore, Eq. 
( [7^ ) remains approximately valid while W drops from being near 1 until it becomes 
sufficiently negative that Y no longer has a small magnitude, though the right hand 
side is no longer correct to first order in the small quantity Y; an expression correct 
to first order in Y for —y/a/Moo 'tiY <^ y/a/Moo, which implies \ W\ <^ 1, is 

^ '''-'l + 2Y)-\ (75) 



Finally, when W is very large and negative, Y can grow from its tiny value and 
gives, for 0.2 < y < -y/a/Moo, 

W M 



1 + sr)-^/^ (76) 



YVl - W Va 



If one only needs an expression for W/ {Y\^l — W) that is approximately correct 
but not necessarily correct to first order in Y, then Eq. (|76D can be taken to apply 
over the whole allowed range where Eq. (|68D is valid, namely —0.2 < Y < Yq 1. 
One can then algebraically solve Eq. (|76D for W = 1 + Sttt'^Pb to get the explicit 
approximate formulas for W = 1 + STrr^Pe and U = 1 — 2m/ r = WY in terms of 
Y = U/W: 



W^^Y(1 + 5Y)-'/' 



2 



Ml 



2a 



-F2(i + 5y)-i4/5^ (77) 



\ a 



4a 



-^F3(1 + 5F)-^^/^ (78) 

One can go on to use Eqs. (|68D-([fO|) to solve for the behavior of r, m, and 0. 
To avoid expressions that are too cumbersome, it is useful to divide up the region 
where Eq. (|68|) is valid into three overlapping regions and give explicit approximate 
results for each region separately. 

In Region 1, where \ — W \ oyY ^fajM^, one gets 

2Moo 

(79) 
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In Region 2, where U <^ 1, which is always true for positive Y <Yo 1 and is 
true also for negative Y if ~Y <^ (M^/aY^^, one gets 



2Moo(l 



4a/Ml + 



m 



1 + 



2^^ 



°°Y^ + 



Moo 




iln. , 

4 Mi 



- sinh ^ — 
2 2Va 



(82) 
(83) 
(84) 



2Moo + 4v^ at F = 0, 
= rj2 ^ Moo + 2^^. 



One can see that r reaches a minimum value of roughly 
where the mass function m attains its maximum value 
Then as one moves further inward, to negative values of Y and W, r increases again, 
and m decreases. 

In Region 3, where —W ^ 1 or —Y ^ ^/a/Moo, one gets 

r f«2Moo(l + 5F)-^/^ 



m 



M. 



iln 
2 



a 

-M^Y 



1 + 5r 



(85) 
(86) 



a 



+ -ln(l + 5r). 

5 



(87) 



One can see that the mass function m crosses zero at F ^ (a/M^)^/^ and thereafter 
becomes negative, approaching negative infinity as Y approaches its lower limit of 
-0.2. 

If we solve Eq. (^) for Y in terms of r, then we can write the approximate 
metric in Region 3 in terms of explicit functions of r: 



5a /2M 



Ml 



/2Mo 



+ 



125a /2Mo 
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Ml V r 
+ r2(d^2^sin'^d(^ 



1 



2M. 



df 

-3 



dr' 



(88) 



One can see that there is a naked singularity, where r goes to +oo and —Qqq goes to 
zero, at a finite proper radial distance. In fact, for r ^ M^o, the proper distance to 
the singularity along a geodesic with constant t, 9, and ip is, if one were to continue 
using this metric in a region where the curvature it indicates is comparable to the 
Planck values, 



'l25a /2Mo 
3 



(89) 
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Therefore, for I -C ^/oi, the metric (P5|) has the approximate form 




1/6 



^^^^ {2M^fr^/\de^ + sm^ 



'■)■ 



(90) 



A shghtly cruder form for the approximate metric of the Boulware vacuum region 
with backreaction, but one which apphes over the whole spacetime, uses (in this 
section only; elsewhere R is the circumference of the shell divided by 27r) 



R 



2M. 



2M. 



(91) 



as the independent radial variable, since this variable, like e^"^ = — (^oo? varies mono- 
tonically with radial distance, from R = 2Moo at the naked singularity to i? = cxo 
at radial infinity. Then the metric takes the approximate form 



1 - 



2M. 



+ 1 



+ R' 



R 

8a 

oo 

1 + 



dt^ 

R 



dR^ 



Qa 



1/3 



{de^ + sin^ 



(92) 



Ml{l-2M^/R)_ 

For a/M^ <^ 1 as we have always been assuming, this metric reduces to very 
nearly the Schwarzschild metric for R — 2M^ S> a/M^, which includes the region 
with Y ^ Yq for which Eq. (|68D is not valid, and it also gives a reasonably good 
approximation to the metric in Regions 1 and 2 where Eq. (|68D is likely to be valid. 
In the part of Region 3 where r ^ Moo, the metric (|92D gives a proper distance 
to the singularity that is roughly VO.9 times the distance in the metric ( pOD at 
the same value of — (^oo; and the circumference of the spheres in the metric (^) is 
roughly (1.2)^/^ the amount given by the metric (pO]), but at least the qualitative 
behavior agrees. Furthermore, the metric (^) is supposed to apply at a distance ^ 
that is less than one Planck length from the naked singularity, and there one would 
expect quantum gravity effects to change the form of the metric or invalidate the 
use of a semiclassical metric altogether. Therefore, I propose that the metric ( P^D 
may be taken as a reasonably good approximation to the metric of an asymptoti- 
cally flat static spherically symmetric vacuum region when the backreaction of the 
Boulware stress-tensor is self-consistently taken into account in a semiclassical ap- 
proximation, and when one avoids the high curvature region where the semiclassical 
approximation is expected to break down. 

We found in Eq. (|55D that for a neutral spherical black hole of area A = AirM'^ 
surrounded by a perfectly reflecting shell at i? — 2M <C M, the entropy is roughly 



32naR 
R-2M 



4 



360(1 - 2M/R)' 



(93) 
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The last term represents the leading term for the reduction of the entropy below 
one-quarter the area. Let us ask how large this term can be for various assumptions 
about the shell. 

First, consider the case that the shell is held up entirely by its own stresses, with 
no external forces (other than gravity) on it. In particular, we shall consider the 
static shell junction conditions (|38D and (p9D, applying the strong energy condition 
to the shell so that its surface stress obeys the inequality 5*1 = 5*1 < —Sq. As we 
shall soon see, it then turns out that f/_ = 1 — 2M_/_R ^ 1 — W = 1 — 2M/ R cannot 
be very small, so the terms involving the pressures inside and outside the shell are 
then negligible. Then the strong energy condition applied to the junction conditions 
([38| ) and (|39D imply that U^U^ > 1/25, and since Eq. (|38|) implies that a shell with 
positive local mass has <U-, we see that 1 — W ^ U- > 1/5, ori?> 2.5M. If 
Eq. (p3D applied for such a large value of 1 — W, it would then give 

-A-S^ ^ , , < §— ^, (94) 

4 360(1 - 2M/R) 72 ^ ^ 

a quite negligible decrease in the entropy, unless somehow nj, + (7/8)raj is very large. 

Next, consider the case that the shell has charge Q, so that its electrostatic 
repulsion holds it up. Since we found above that the stresses within the surface 
of the shell are quite ineffectual in holding up the shell at R — 2M -C M, let us 
drop them from the junctions equations but add the tension of the electromagnetic 
field outside the shell and assume that that tension is much greater than the radial 
pressures (or tensions) of the quantum fields. Then the junction conditions (^) and 
( p9| ) become 

^ = U'J' - UT (95) 
R 



= 8F = ^ + U:'^' - (1 - %)U+'^\ (96) 



and 

= HF = 

R ' R^ 

1 /2 

Now for fixed charge-to-mass ratio Q/yU, if we let 7 = {fi/R)/U_ < 1, Eq. ( pS]) 

1 /2 

imphes that Uj_ = (1 — 7)(/x/i?), which when inserted into Eq. ( p6D gives 

7^ = 1 - 7 + 7(Q//^)' < (Q//^)'. (97) 



1 - 2M/R U. 

If we take the charge-to-mass ratio of an electron, we get {Q/^Y ~ 4.17 x 10^^. 
If we then suppose that somehow a shell of electrons reflects electromagnetic (but 
not other) radiation and thereby manages to keep the electromagnetic field in its 
Boulware vacuum state outside the shell (rather than in the Hartle-Hawking thermal 
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state that exists within the shell), then Ub = 2 (from the two helicities of photons) 
and Uf = 0, so one gets 




1 



,40 



(98) 



180(1 - 2M/R) 



< 2.31 X 10 



Of course, there are severe problems in attaining anything near this limit. First, 
electrons in a shell around a black hole, even if in static equilibrium as I have calcu- 
lated they can be, will not be in stable equilibrium, and some unknown mechanism 
would have to be invoked to keep the shell in place. Second, without specifying 
how the electrons are to be kept in place, it is hard to say how they will respond to 
the black hole thermal radiation impinging upon them from below. However, it is 
interesting that the upper limit given by Eq. (^8]) for the reduction in the entropy 
below one-quarter the (neutral) black-hole area, from a shell held up by electrostatic 
forces, is so large (because the charge-to- mass ratio of an electron is so large). 

For a somewhat more nearly realistic example of a shell around a black hole, 
consider a thin aluminum foil that is charged so that, like the shell of pure electrons, 
the electrostatic forces balance the gravitational forces. In this case there will be 
limitations from the mass density p of the foil, the minimum practical thickness r 
of the foil and the maximum charge per surface area, a, that it can hold. Again 
the possible tensile forces within the foil itself are negligible in comparison with the 
electrostatic forces and hence will be ignored. (It turns out that the strong energy 
condition also implies that they are not nearly sufficient to stabilize the shell against 
radial perturbations, which are unstable because the local gravitational forces go up 
rapidly as the shell is brought closer to the black hole horizon, whereas the electro- 
static repulsion forces depend only on the circumference of the shell, which changes 
only slowly as the shell is moved inward or outward near the horizon. However, 
just as modern jet fighter planes fiy under the control of computer servomechanisms 
while being deliberately constructed to be unstable, and thus rapidly maneuverable, 
I shall suppose that some unspecified servomechanism can be used to keep the shell 
in place. I shall leave to the reader the engineering problem of constructing such a 
servomechanism and just tell how to balance the forces in the unstable equilibrium.) 

Let me now give some parameters associated with the aluminum foil, using both 
conventional units, atomic units, and Planck units (always the case in this paper 
when no units are explicitly given). For this discussion I shall use the charge of the 
positron as 




(99) 



the mass of the positron or electron as 




(100) 
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the mass of the proton as 



mp = mp/^Tic/G ^ 7.684 x lO'^o, (101) 

the Rydberg energy as 

En = ^me^ « 13.6057 el/ ^ 1.114 x 10"^^ (102) 

the Bohr radius as 

qb = ^ 5.2917721 x lO^^cm ^ 3.2755 x 10^^, (103) 
me"' 



the tropical year as 
and the solar mass as 



1 yr ^ 5.854 x 10^", (104) 



Mq ^ 9.137 X 10^^ ^ 0.5395/mJ. (105) 
The density of solid aluminum is then 

p = NpTUp/al ^ 2.70 g/cm^ ^ 5.233 x 10'^^, (106) 

with 

Np ^ 0.2391. (107) 

The aluminum will be so cold it will be superconducting, in which case it has a 
Meissner magnetic penetration depth of about 50 nm As we shall see, the local 
Hawking temperature will be far below the superconducting gap energy, so the foil 
will be almost completely reflecting to the thermal electromagnetic radiation if its 
thickness is several times the magnetic penetration depth. To be very conservative, 
I shall take the thickness of the foil to be about 100 times the magnetic penetration 
depth, 

r = NrttB = 50 microns = 0.0005 cm ^ 3.094 x 10^^ (108) 

with 

AT^f^ 94486. (109) 

I shall choose the electric surface charge density a so that if it were an excess of 
electrons, the probability for one to tunnel off, say P, is very small, say exp (—100) 
in some suitable atomic time unit. Since P ~ exp (—2/) with tunneling amplitude 
/, I shall choose / = (1/2) In (1/P) = 50. Now the work function for polycrystalline 
aluminum is 

Vo = vqEr ^ A.28 eV, (110) 
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with 

t;o^ 0.315. (Ill) 

If one takes the potential energy of an electron, relative to that at the Fermi surface, 
to be Vq — eEx at distance x from the surface, where E = Aira is the magnitude of 
the charge density (dropping the minus sign that a would have in actuality, because 
of Ben Franklin's inconvenient sign convention, if there were an excess of electrons 
on the surface), then the WKB amphtude for the electron to tunnel through the 
classically forbidden region < x < xq — Vo/{eE) is 

fxo , \ SmVJf 

I = / \/2mVdx = ^ ^ ^ . (112) 
Jo 6eE 



This then gives 



•2 'S 2^3/2 2^3/2 

m e m e Vn m e Vq 

a = 



67rln(l/P) GOOtt 
3.34 X 10^^ e/cm^ ^ 0.00893 e/a| 7.46 x 10"^^ (113) 



with 



= 67rln {l/P)vo^^^ = 6007rvo"^^^ 10 700. (114) 

Expressed in terms of the unit area n~^^^ formed from the aluminum atomic number 
density n, one needs about one excess electron per 460 of these unit areas to give 
this surface charge density a, so this does not seem excessive. 

From these parameters, one gets that the charge-to-mass ratio of the aluminum 
foil is 

Q a e/m„ 1.112x10^^ 

— = — = — '-^_r^ — fti4.61xl0^ (115) 

II pr N^NpNr 2.41 x lO^ ^ ' 

This is that of the pure electron shell, divided by NaNpNrimp/m) ^ 4.43 x 10^^ 

Then, by the same analysis used above for the pure electron shell, one finds that 
if one takes ji/R— U]!"^ — n/Q, one gets 

^ 2.12 X 10l^ (116) 



1 - 2M/R \ II 

and hence the reduction of the entropy from excluding the thermal photons from 
above the shell is 



4 180(1 - 2M/R) im{N^NpNrf 

1.24 X 10^6 



1.05 X 10^9 



1.18 X 10^^ (117) 
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As we shall see, this is small in comparison with the total entropy of the black hole, 
but the reduction means that the number of states is fewer by a factor of about 



AS 1 n51 000 000 000 000 000 



(118) 



which is quite a large factor. 

One can further derive properties of the black hole around which this aluminum 
foil is placed to minimize 1 — 2M/R and hence maximize AS*. The radius of the 
shell (which is very nearly that of the hole) is 

R = ^ 1.07 X 10^^ ^ 182 light years, (119) 
giving a circumference of about 1150 light years. If the redshift factor to infinity were 



1 — 2M/ R, the time as seen at infinity for a photon to circumnavigate the shell 
would be about 5.3 x 10^^ years, hundreds of times longer than the current age of 

1 /2 

the universe. (In actuality, if one does take the limit of setting fi/R = UJ = fi/Q, 
one finds that as seen from the outside, the shell is at what would be the horizon 
of an extreme Reissner- Nordstrom black hole, so there would be an infinite redshift 
factor to infinity. But if one set fi/R to be, say, half as large, the previous quantities 
would be shifted by merely factors of two or so, and then the time as seen at infinity 
for a photon to circumnavigate the foil shell would be of the order of hundreds of 
times the present age of the universe.) 
The mass of the black hole is then 

M ^-R = ^ 5.34 X 10^2 ^ 5 84 ^ lO^^M©, (120) 
2 Ana 

of the order of mass of a supercluster of galaxies. This mass then gives a total 
entropy for the black hole of 

S ^ AnM^ ^ 3.58 x 10^°^ (121) 

which is about 3.04 x 10*^ times the reduction AS" in the entropy calculated above. 
Therefore, as already mentioned, the relative reduction of the entropy is negligible 
in this case, but it does reduce the total number of states (far, far more than a 
googolplex in this example) by the huge factor given by Eq. ( |118| ). 
The proper radial distance of the shell to the horizon is 



D ^ 2M^1 - 2M/R ^ 2.50 light seconds. (122) 
The local acceleration of gravity as seen from just inside the shell is 

^ ~ _ ~ 1.20 X 10* m/s^ ^ 1.22 x 10^ g^, (123) 
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over twelve million times the acceleration of gravity ^f® at the surface of the earth. 
The local Hawking temperature as seen at the inner edge of the shell is 



Tiocai ~ ~ 3.44 X 10-^5 ^ 4.87 x 10-^=^ K, (124) 



which is indeed far below the critical superconducting temperature of 1.175 K | pi][ , 
so the shell should stay superconducting and indeed reflect almost all of the electro- 
magnetic radiation emitted by the black hole. 

Thus we have seen that by placing a reflecting shell around a black hole, we can 
make the entropy have a value that is below one-quarter its area. If we are allowed 
an idealized perfectly reflecting shell that can be placed within roughly one Planck 
length of the horizon, then this entropy reduction can be of the same order as the 
area of the hole. For a more realistic shell, such as a superconducting aluminum 
foil, the entropy reduction can only be a tiny fraction of the area, but it still can be 
huge in absolute units, markedly reducing the number of black hole states from what 
would be erroneously estimated by exponentiating one quarter the horizon area. 
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not immediately come to mind, Valeri Frolov, Frank Hegmann, Akio Hosoya, Satoshi 
Iso, Sang-Pyo Kim, Frank Marsiglio, Sharon Morsink, Shinji Mukohyama, Jonathan 
Oppenheim, L. Sriramkumar,and Andrei Zelnikov. Part of this work was done at 
the Tokyo Institute of Technology under the hospitality of Akio Hosoya, and part 
was done at the Haiti Children's Home of Mirebalais, Haiti, under the hospitality 
of Patricia and Melinda Smith while adopting six-year-old Marie Patricia Grace to 
become our fifth child. This research was supported in part by the Natural Sciences 
and Engineering Research Council of Canada. 
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